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Abstract 

We show that Brownian motion is spatially not symmetric for mesoscopic particles embedded in a fluid if 
the particle is not in thermal equilibrium and its shape is not spherical. In view of applications on molecular 
motors in biological cells, we sustain non-equilibrium by stopping a non-spherical particle at periodic sites 
along a tilament. Molecular dynamics simulations in a Lennard- Jones fluid demonstrate that directed motion 
is possible without a ratchet potential or temperature gradients if the asymmetric non-equilibrium relaxation 
process is hindered by external stopping. Analytic calculations in the ideal gas limit show that motion even 
against a fluid drift is possible and that the direction of motion can be controlled by the shape of the particle, 
which is completely characterized by tensorial Minkowski functional. 
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Mesoscopic particles dissolved in a fluid are expected to perform symmetric thermal fluctu- 
ations around a mean position . The shape of the particles is irrelevant for this so-called 
Brownian motion - as long as the particle is in thermal equilibrium with the fluid. Consequently, 
a net transport in a preferred direction is not possible without applying an external force which 
breaks the spatial symmetry. In the past, several models for Brownian motors have been proposed 
based on different methods of rectifying thermal noise; for instance, by an asymmetric external 
potential which is switched on and off periodically (for a review see Ref. U). Another concept 
is based on two heatbaths at different temperatures such as the Feynman ratchet [4] and its sim- 
plification by Van den Broeck ISO. It was shown that directed Brownian motion can be achieved 
as long as the two reservoirs have different temperatures and spatial symmetry is broken in some 
way. This letter demonstrates that a Brownian motor can be built even in a single heatbath without 
violating the second law of thermodynamics. Hence, it is a possible theoretical model for molec- 
ular motors in biological cells. We show that a steady state motion of non-spherical particles is 
obtained if their relaxation towards equilibrium is prohibited by periodical stopping. 

To illustrate the concept we restrict the system to two dimensions and follow Ref. [pO. Since 
the spatial symmetry breaking is a prerequisite for net transport, we consider an asymmetric 
motor K of mass M in a single heatbath at temperature T. The motion of the motor is re- 
stricted to a one-dimensional track, here the x-axis, so that the motor velocity can be written 



as V = {V^JkBT/M , 0) with the normalized velocity V in x-direction. Furthermore, elastic in- 
teractions are assumed between the motor K and the surrounding fluid particles. In contrast to 
Refs. itzilsl] there is no dissipation involved in this collision. However, due to the confinement only 
the axial momentum component is conserved. In thermal equilibrium, directed motion in a single 
heatbath is prohibited by the second law of thermodynamics. The motor K behaves like a Brow- 
nian particle with a Maxwellian velocity distribution with mean (y) = and variance {V'^) = 1; 
independent of its particular shape. Hence, in order to get directed motion, we have to sustain a 
non-equilibrium state. 

Our model is motivated by a molecular motor in a cell, e.g. kinesin, that binds to a filament 
after every step [l6D. Similarly, our particle is stopped at periodically spaced binding sites along the 
track, i.e. its kinetic energy is set to zero whenever it reaches a binding site. Thermodynamically, 
of course, the stopping requires a greater amount of work than the decrease in entropy; here ^/c^T 
from the reduction of the Maxwellian distribution for V to the velocity distribution PiV) = 5{V) 
for the stopped particle. After binding the motor is released again and starts to relax towards 
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thermal equilibrium. By means of collisions between fluid particles and the motor, kinetic energy 
is transfered from the fluid to the motor until the motor has the same effective temperature as the 
surrounding bath. 

In the following, we show that this non-equilibrium relaxation process is asymmetric due to 
the asymmetry of the particle. This asymmetric relaxation process can be used to rectify thermal 
fluctuations. To test this hypothesis, we performed 2D molecular dynamics simulations of an 
asymmetric particle in a two-dimensional fluid. The interaction between the fluid particles is 
described by the Lennard- Jones (L J) potential W(r) = 48 (^(7)^^ — (7)^^ where r is the distance 
between two particles with diameter a and £ is the depth of the potential well. In the following, 
all lengths are given in units of a, energies in units of £ and times are expressed in units of to = 



(7^/m/£ where m denotes the mass of the fluid particles. The simulation box with edge length 
lOOcr contains n = 800 fluid particles yielding a particle density g ^ O.OScr^^. The temperature is 
setto T = 3£ /ks and the overall simulation time was 2000)f:o- With these parameters the LJ fluid 
is in the gas phase without any long range correlations. Periodic boundary conditions are used and 
also a cut-off radius of 2.5a in U{r) to limit the range of the potential. 

As illustrated in Fig. [U the motor is built of stiffly-linked fluid particles which interact with 
the surrounding fluid also by a Lennard- Jones potential. Consequently, the motor mass M = Nm 
scales with the mass m of the fluid particles. The motion of the motor is confined to a one- 
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FIG. 1: An asymmetrically shaped motor (here a triangle) is built from fluid particles and placed in a 
Lennard- Jones fluid. Its motion is restricted to a one-dimensional track with periodically spaced binding 
sites along the x-axis. If the motor's center of mass crosses a binding site, the velocity of the motor is set 
to zero. To test the importance of the motor's shape for the momentum transfer during collisions with fluid 
particles we orient the triangle in two different ways. 
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FIG. 2: Histograms of end positions at time t = 2000to of a triangular motor where the distance between 
stopping sites is 6s (see Fig. [T]). The asymmetrically oriented motor moves preferable to positive X- values 
whereas the symmetric triangle diffuses around X = 0. A discrete spacing is visible since the motors 
accumulate around the stopping sites. 400 realisations are used with integration step 6t = O.OOlto and 
2000000 steps. With decreasing distances 6s between stopping sites the distribution becomes narrower, as 
expected for a hindered diffusion process. 

dimensional track with periodically spaced stopping sites. Every time the motor's center of mass 
reaches one of those sites, the velocity of the motor is set to zero and the corresponding kinetic 
energy is re-transferred to the fluid. 

Consider a motor with the shape of an isosceles triangle of height lOa and basis length 20a 
which corresponds to an apex angle of 6^ = 0.25268, = 53 and mass M = 53m. The motor is 
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placed in the beginning (t = 0) at the position X = 0. Figs. |2] and |3] show the simulation results 
for the motion of the motor. Note, that for the chosen parameters the typical time between two 
successive stops of the motor is smaller than the relaxation time into equilibrium. In case of an 
isosceles triangle, the motor can be oriented symmetrically or asymmetrically with respect to its 
direction of motion (Fig.[T]), to test the importance of the shape. Fig.[2]presents the distribution over 
multiple runs of the end positions, i.e. the positions of the motors after time 2000to, for different 
stopping distances 5s. If the triangle is placed with the apex onto the track, i.e. asymmetric with 
respect to the direction of motion, the maximum of the distribution of the end positions is shifted to 
positive values x > 0. Hence, on average, the triangle travels a certain distance towards its basis. 
This is in accordance with the results shown in Fig. [3] where a net speed of about (v) = 0.0 11 1^ in 
the direction of the triangle basis can be observed for asymmetric orientation. In contrast, aligning 
the triangle basis parallel to the x-axis conserves spatial symmetry, i.e. mirror- symmetry with 
respect to the y-axis, and therefore causes the mean position to fluctuate around zero without any 
net motion (see Fig.|3]). The end positions of the motors are in this case symmetrically distributed 
around the origin (see Fig. |2l). This shows two things: first, the proposed model demonstrates the 
feasibility of rectifying Brownian motion in a single heatbath by stopping the particle at periodic 
sites. Second, it is sufficient to break the spatial symmetry by the shape of the particle, so that an 
asymmetric external potential is not necessary. 

To get a better insight into the microscopic mechanism of the driving force, we calculate ana- 
lytically the stochastic dynamics of this non-equilibrium relaxation process for a simplified model. 
The interaction between the fluid and the motor is no longer a Lennard- Jones potential but an ideal 
hard sphere interaction with infinitesimal interaction time. The total energy is conserved as well 
as the momentum in x-direction and tangential to the particle's surface. Details of this model can 
be found in Ref. dSfl. 

Furthermore, we assume an infinite heatbath, so that the temperature is constant and the velocity 
distribution of the fluid is a Maxwell distribution = ^j;^^ exp (^—2^^{v — . To 

study motion against a drift we introduce an average drift velocity u of the fluid. In the beginning, 
we assume u = 0. In order to exclude multiple collisions, we use the ideal gas limit of the fluid 
and confine the motor shape to convex objects K so that the stochastic motion is a Markov process. 
This allows us to describe the time dependence of the velocity distribution P(y, t) of the motor by 
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FIG. 3: The average position (400 runs) of an asymmetrically oriented, isosceles triangle of height lOcr 
and basis length 20a moves with a constant velocity in the direction opposite to its sharp edge (i.e., in the 
direction of its basis). Thus, directed motion is possible when the triangle is oriented asymmetrically with 
respect to the direction of motion. Otherwise the motor performs symmetric fluctuations around its starting 
position. Inset: trajectories of individual motors. Note, that the velocity decreases for decreasing distances 
5s between stopping sites if the typical stopping time is smaller than the relaxation time of the motor. 

a Master equation with a transition rate [|50 



W{V'\V) = I dO 

'dK 



X 



dv p ^q{v) iy — v) ■ 



(1) 



(y = ^JhsT /M V) which equals the number of collision processes from V to V per unit time 
using the dimensionless time t = Se^J ^^i. Here, ey = (cos 0, sin 0) is the tangential and 



= (sin 0, — cos 0) the normal vector in the collision point r G dK on the boundary dK of 
the motor and is the polar angle. The momentum transfer of a single collision is then given by 



Ap 
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\ — v) ■ e± (see Ref. [5]). Because the Master equation cannot be solved 
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exactly, we apply a Kramers-Moyal expansion yielding 

dP{V, t) 



dt 



^l^.[-w) iMV,t)P{v.t)] (2) 

n=l ^ ' 



withthe coefficients A„(\/,t) = / dV'{V'-V)''W{V\V'). Different to Ref. iSD, we emphasise two 
important points: (i) the dependence on the shape of the motor can be fully described by tensorial 
Minkowski functionals My'^"^ {K) of the motor K; (ii) the non-stationary solution of Eq. ^ is the 
essential ingredient to rectify Brownian motion if relaxation is prohibited by stopping sites (for 
details see Ref. M)- 

Tensorial Minkowski functionals are defined as surface integrals in d dimensions lioll 



s— times 



Mi^'^'Xk) = dO v0...0r0e±0...0e± (3) 

JdK 

over tensor products of position r and normal vector on the surface dK of motor K. Us- 
ing (e^ ■ e^)" = sin" 0, the x . . . x-components of the first tensorial Minkowski functional in 

/ , \ n—times 

two dimensions can be written as (Mr"^(ir))^^^ = Jgj^dO sin" (h. For convex bodies 
K one finds mI)^'^\k) = so that the first two jump moments read Ai{V) = r^^iV + 
4Kiax(V"^ - 1)) and A2(y) = 2r"^(l + 12V^!,J^). Here, we introduced the relaxation rate 
r"^ = y/p (m[^''^\k) j /2 and the maximum velocity 

V /XX 

K.ax = - a/tT ■ (4) 



4 V8 (Mr)(K)) 

V / : 



The shape dependence of the stochastic dynamics enters only via these two parameters which 
can be expressed in terms of Minkowski tensors. These functionals are known in many cases to 

o 

allow for a complete description of the shape-dependence of physical properties [|10|1 . The other 
Minkowski tensors in Eq. ([3]) play a role when the collision rules in Eq. ([T]) are modified or the 
drift is non-zero, for instance 

From Eq. one immediately obtains the time dependence of the average velocity dt{V) = 
{Ai{V)) and the average squared velocity dt{V^) = 2{VAi{V)) + {A2iV)). The equations are 
coupled to higher moments of V because of the dependence of AniV) on V. However, assuming 
the mass m of the fluid particles to be small against the motor mass M, an expansion in e = a/p' 
leads to a decoupling of the differential equations from higher moments when terms of order 0{e^) 
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{V)+ {V')-l\V, 



or higher are neglected 

d{V) _ _1 
dt ~ ~T 

'-^ = l{l-{V^))+0{e^) . (5) 

Hence, the solutions give readily the mean velocity {V){t) = 4Vniax (e~*/^ — e~^*/^), (V^)(t) = 
1 - and the average of the travelled distance = dt' V{t') = 2V^^^t (l - e-*/^)l 

The time r characterizes the increase of the kinetic energy ^{V'^) towards its equilibrium value 
^ksT, which is the stationary solution of Eq. ([5]). The relaxation process depends on the shape 
of the motor only via the maximum of the averaged velocity V^ax^ so that the shape depen- 
dence is fully described by the Minkowski functionals of the particle K. For instance, the 
Minkowski functionals of a triangular motor with the top pointing in negative x-direction reads 
(Mj°'''(ir)),, = L(l + sin^)and (M^ '''(iT)),,, = Lcos^^, so that K.ax = (l-sin0) 
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and = Ly/p{l + sin 9)/2 depend on the apex angle 9 of the triangle (see Fig.[T]). 

The shape dependence of this non-equilibrium Brownian motion can be perfectly demonstrated 
with the 'capped-triangular' motor shown in the inset of Fig. Sfb), whose Minkowski functionals 
are (M^ '^^(K))^^ = ^ (f + sin^) and (Aff '^^(/sT))^^^ = ^ (1 " sin^ 6). As shown in Fig.Htb), 
the relaxation time r decreases with increasing opening angle 9 whereas the maximum average 
velocity K^ax increases with 9 from negative values to positive ones. Interestingly, one can change 
the direction of motion just by adjusting the shape of the motor. For small 9, i.e. long triangle legs 
relative to the arc length of the semi-circle, the motor moves towards the semi-circle, whereas it 
travels towards the triangle if 9 is close to |. Note, that there is exactly one value for 9 for which 
the motor does not move on average (Kiax = vanishes) though the motor definitely breaks spatial 
symmetry. 

The relaxation process of the triangular motor is shown in Fig. SI The average velocity first 
increases until a maximum value l^^iax is reached and than decays exponentially towards the equi- 
librium with vanishing mean velocity. Note, that in the beginning the motor travels an average 
distance X^ax = 2KnaxT before it fluctuates around the mean position X^ax in the long time limit. 
Shortly after releasing, the motion of the motor has a preferred direction. Thus repeated stopping 
yields directed motion if it occurs before the system has relaxed completely. 

A simple analogy of the rectification mechanism is a piston in a cylinder separating two fluids 
in thermal equilibrium at the same temperature. If the piston is arrested or has an infinite mass, 
one finds in thermodynamic textbooks that the averaged momentum transfer on one side per unit 
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FIG. 4: (a) A triangular motor (apex angle = vr/G, base length = 1, = m/M = 0.05) can swim 
even against a fluid drift velocity U = -y/ M/ {UbT) u during its relaxation in thermal equilibrium, (b) The 
direction of motion can even be inverted solely by changing the shape of the motor. The maximum velocity 
14nax and the relaxation time r of a 'capped-triangular' motor (sketched in the inset) can be regulated by 
the opening angle 9 which determines the Minkowski functionals in Eqs. ([3]l and 

area of its cylindrical cross section equals the pressure of the fluid on this side. Since the fluids 
on both sides are equilibrated at the same temperature there is no net force on the piston. The 
same is true if the piston can move freely and is in thermal equilibrium with the fluids on both 
sides. However, a piston with finite mass M and not in equilibrium experiences kinematic effects 
due to momentum conservation when momentum is transfered by collisions, so that the net force 
depends on the shape of the piston. 

The motor is also able to move against small fluid drifts. Repeating the analytic calculations 
with tt 7^ and confining the drift to the s-direction we find the mean velocity shown in Fig. |4ta) 
for a triangular motor with its apex pointing in negative x-direction. In the long time limit the 
motor equilibrates always to the drift velocity independent of sign and strength of the drift. For 
small drifts in the preferred direction - here the positive x-direction - the motor accelerates to a 
maximum average velocity V^ax above the equilibrium value before it relaxes to the drift velocity 
u. When applying a small drift in opposite direction, it can be seen that the motor starts to move 
in positive direction but changes its direction of motion after a certain time. When stopping it 
repeatedly before the direction is changed, an average motion against a fluid drift is achieved. 



9 



1.6 



0.8 



0.2 



0.4 



0.6 



1.4 



1.2 



f 



1 




15 



5 



10 



20 



t/t 







FIG. 5: Relaxation of a motor to thermal equilibrium: the molecular dynamics simulations of the LJ-fluid 
(symbols; averages over 400000 runs) can be well described by the analytic results in the ideal gas limit 
(hues) if the data are normahzed by the maximum velocity Knax and the relaxation time. 

Fig.|5]compares our analytic results with the numerical data from the molecular dynamics simu- 
lations. Evidently, a real fluid with Lennard- Jones interactions cannot be in quantitative agreement 
with calculations based on an ideal gas. Nevertheless, we find the same functional form for the 
non-equilibrium relaxation process given by the solution of Eq. dS]), which confirms our theoreti- 
cal understanding of this rectification mechanism. We showed that thermal fluctuations in a single 
heatbath can be used to get directed motion if (i) the motor is sustained in a non-equilibrium state 
by stopping it at periodically spaced sites on a one-dimensional track (or alternatively at periodic 
time intervals); (ii) the spatial symmetry is broken by an asymmetric shape of the motor which 
makes the non- stationary relaxation process also asymmetric. The shape-dependence of the non- 
equilibrium relaxation is used to get directed motion even against small fluid drifts. By changing 
the shape of the motor, it is possible to change the direction of motion and to control the average 
velocity of the motor. 
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